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O '■ Abstract 

O 

o, 

CN , We derive analytic expressions for three flavor neutrino oscillations in the 

^ I presence of matter in the plane wave approximation using the Cayley- 

|JZ ' Hamilton formalism. Especially, we calculate the time evolution operator in 

(<— ^ ■ both flavor and mass bases. Furthermore, we find the transition probabilities, 

cn . matter mass squared differences, and matter mixing angles all expressed in 

.^ , terms of the vacuum mass squared differences, the vacuum mixing angles, and 

^ I the matter density. The conditions for resonance in the presence of matter 



are also studied in some examples. 



I. INTRODUCTION 



o 
o\ 

■^ I There are at present essentially four different types of experiments looking for neutrino 

p ■ oscillations [|I|: solar neutrino experiments, atmospheric neutrino experiments, accelerator 

A . neutrino experiments, and reactor neutrino experiments. Among the accelerator and reactor 

D ! experiments, the accelerator long baseline experiments [Q] are just beginning to become 

'^ ' operative. Also the Sudbury Neutrino Observatory (SNO) @] has started to take data 

on solar neutrinos. An accumulating amount of data on neutrino oscillations is therefore 
/\^ • becoming accessible. Yet, at present there is no general agreement on the mixing angles or 

j^ . the mass squared differences of the neutrinos. 

In a previous paper 0] , we described a global analysis of neutrino oscillation data in a 
three flavor neutrino scenario. In that analysis, we deliberately considered mass ranges for 
the solution that avoided the region where the solar and atmospheric neutrinos could be 
affected by the so-called Mikheyev-Smirnov-Wolfenstein (MSW) effect |^. 

Here we would like to report on the second stage of a global analysis of all the neutrino 
oscillation data in the three flavor neutrino model and present analytic expressions for the 
oscillation probabilities in the presence of matter expressed in the Cabibbo-Kobayashi- 
Maskawa (CKM) mixing matrix elements (or [vacuum] mixing angles) and the neutrino 
energies or masses, i.e., incorporating the MSW effect. 
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We will assume as before that CP nonconservation is negligible at the present level of 
experimental accuracy [^]. Thus, the CKM matrix for the neutrinos is real. 

Previous work on models for three flavor neutrino oscillations in matter includes works 
by Barger et al. 0, Kim and Sze [Q, and Zaglauer and Schwarzer 0. Our method is 
different from all these approaches and their parametrizations also differ slightly from ours. 
In particular, we calculate the time evolution operator and do not introduce the auxiliary 
(effective) matter mixing angles. 

Approximative solutions for three flavor neutrino oscillations in matter have been pre- 
sented by Kuo and Pantaleone |TD] and Joshipura and Murthy [^. Approximative treat- 



ments have also been done by Toshev and Petcov |T2|. D'Olivo and Oteo have made contri 



butions by using an approximative Magnus expansion for the time evolution operator fl^ 



and Aquino, Bellandi, and Guzzo have found the auxiliary matter mixing angles in terms 



of the vacuum mixing angles [0. Extensive numerical investigations of matter enhanced 



three flavor oscillations have been made by Fogli et al. [IT3 



Below we derive the expressions for the transition probabilities in the case of three fla- 
vor neutrino oscillations when the neutrinos pass through matter (the MSW effect) using 
Cayley-Hamilton's theorem. This allows us to give rather explicit analytic formulas in terms 
of elementary functions for the quite involved expressions for the transition probabilities. 
We derive the time evolution operator for a constant matter (electron) density. The most 
straightforward way to use this result is to approximate any density profile with step func- 
tions and to obtain the total time evolution as a product of the evolution in each step. This 
can easily be handled numerically |jT6[ . 



The three auxiliary angles of rotation from the matter mass eigenstate basis to the flavor 
state basis can be identified from these expressions and we show their relation to the ordinary 
formalism. 

The application of this formalism to the analysis of the neutrino oscillation data will be 
presented in a forthcoming publication. Here we will only illustrate the qualitative behavior 
of the transition probabilities and the auxiliary matter mixing angles as functions of the 
matter density for different choices of the vacuum mixing angles and the mass squared 
differences. 

It is possible to adapt the formalism presented here to the case of a scenario with an 
arbitrary number of neutrinos, but the eigenvalues have to be solved for numerically for the 
case of five or more neutrino flavors. 

The outline of our paper is as follows. We first go through some formalism in Sec. |I| 
and then set out to find the solution for the time evolution operator and the transition 
probabilities in the presence of matter. In Sec. Q, the main result of our analysis is given 
by the time evolution operator for the neutrinos when passing through matter with constant 
matter density expressed as a finite sum of elementary functions in the matrix elements of 
the Hamiltonian, and in Sec. |V|, the transition probabilities in presence of matter are given. 
In Sec. 0, we determine the auxiliary matter mixing angles. Finally, in Sec. 0, we present 
a discussion of our main results and the illustrations. 



II. FORMALISM 

Neutrinos are produced in flavor "eigenstates" \va) (« = e, /i, r) created by the interaction 
of tlie weak gauge bosons witli tlie cliarged leptons. Between tlie source, tlie production point 
of tlie neutrinos, and the detector, the neutrinos evolve as mass eigenstates |z/(j) (a = 1, 2, 3), 
i.e., as states with definite mass. At the detector, the neutrino flavors can again be identified 
by charge current interactions. It is also possible to detect the neutrinos by weak neutral 
current interactions. In this case, however, the weak neutral currents do not distinguish the 
flavors on the detector side, and the total neutrino flux will be measured. 

Let the flavor state basis and the mass eigenstate basis be given by Tij = {\i^a)}a=e,^i,T 
and Tim = {Wa)}l=i, respectively. Then the flavor states l^a) G 'Hf can be obtained as 
superpositions of the mass eigenstates |z/a) G TCm, or vice versa. Observe that the bases TCf 
and Tim are just two different representations of the same Hilbert space Ti. 

In the present analysis, we will use the plane wave approximation to describe neutrino 
oscillations. In this approximation, a neutrino field u^ with flavor a is a linear combination 



of neutrino mass fields z/^ such that [0 



/„ = ^ UaaJ^a, a = 6, fJ., T, (l) 



a=l 



where the UaaS are entries in a 3 x 3 unitary matrix U. Taking the Hermitian conjugate of 
Eq. (D yields 



3 
a=l 



XI ^""'^a' a = e,/i,r. (2) 

a=l 

Acting with Eq. (|) on the vacuum state |0) gives 

3 

'a) = ^KaWa), a = e,fi,T, (3) 



a=l 



since the states are defined as |z/q) = z/||0), where a = e,fi,T, and |z/a) = z/^jO), where 
a = 1,2,3. In what follows, we will use the short-hand notations |a) = |z/q,) and \a) = {ua) 
for the flavor states and the mass eigenstates, respectively. An arbitrary neutrino state 
^ G 7i can, of course, be expressed in both the flavor and mass bases as 

3 3 / \ 3 

i^= Yl V^"l«)= 5Z V'a5Zf/:al«)=$^( Yl ^aU:A\a)=J2^a\a), (4) 

a=e,fi,T ct=e,fjL,r a=l a=l \a=e,^,T / a=l 

where ipa, « = e,fi,T, and ipa, a = 1,2,3, are the components of the state ip in the flavor 
state basis and the mass eigenstate basis, respectively, and they are related to each other by 

^a= Yl ^->-' « = 1.2,3. (5) 



In matrix form, one has 

^/ = U^Jrn, (6) 

where 

^/ = (^a) = I i>i, \ eHf and l/jm = (^a) = ^"2 e Km- 

Here the subscripts / and m denote the flavor state basis and the mass eigenstate basis, 
respectively. 

The unitary matrix U is given by 

U = (U^a) = U^l U,2 U,3 . (7) 

V Url Ur2 UrS J 

A convenient parametrization for U = U{6i, 62, O3) is given by [^ 

(C2C3 S3C2 S2 \ 

—S3C1 — S1S2C3 C1C3 — S'iS'25'3 S1C2 I , (8) 

O1O3 — O2C-1C-3 — O1C-3 — O2O3O1 0102 / 

where Si = sin^j and Ci = cos6i for i = 1,2, 3. This is the so-called standard representation 
of the CKM mixing matrix. The quantities 6i, where i = 1,2, 3, are the so-called (vacuum) 
mixing angles. We have here put the CP phase equal to zero in the CKM matrix. This 
means that ?7*„ = Uaa for a = e, fi,T and a = 1,2, 3. 

In the mass eigenstate basis Tim, the Hamiltonian J€^ for the propagation of the neutrinos 
in vacuum is diagonal and given by 

f El \ 
Hm= \ E2 ], (9) 

V eJ 

where Ea = \/'m\ + p^, a = 1,2,3, are the energies of the neutrino mass eigenstates \a), 
a = 1,2,3 with masses rria, a = 1,2,3. Note that we assume the momentum p to be the 
same for all mass eigenstates. 

When neutrinos propagate in matter, there is an additional term coming from the pres- 
ence of electrons in the matter . This term is diagonal in the flavor state basis Tif and is 
given by 



Vf = A\ , (10) 




where A = ±\^GpNe is the matter density. Here Gp is the Fermi weak coupling constant 
and Ne is the electron density. The sign depends on whether we deal with neutrinos (-I-) or 
antineutrinos (— ). We will assume that the electron density N^ is constant throughout the 



matter in which the neutrinos are propagating. In the mass representation, this piece of the 
Hamiltonian is Vm = U~^VfU, where again U is the CKM matrix. 

We are looking for the unitary transformation that leads from the initial state ipf{0) in 
flavor basis at time t = of production of the neutrino, to the state of the same neutrino 
ipf{t) at the detector at time t. This transformation is given by the operator f//(t) = Uf{t, 0), 
where Uf{t2,ti) is the time evolution operator from time ti to time ^2 in flavor basis. This 
operator can be formally written as Uf{t) = e"*^-^*. When the neutrinos are propagating 
through vacuum, the Hamiltonian in flavor basis is Hj = UHmU~^. In this case, the 
exponentiation of Hf can be effectuated easily: Uf{t) = e~'^^f^ = f/e~*^'"*f/~^, and the 
result can be expressed in closed form. In the case when the neutrinos propagate through 
matter, the Hamiltonian is not diagonal in either the mass eigenstate or the flavor state 
basis, and we have to calculate the operator Uf{t). This is the task we now set out for. 

III. CALCULATING THE TIME EVOLUTION OPERATOR 

The Schrodinger equation in the mass eigenstate basis is 

^irJ^rn{t) = ^J^m^.mit), (11) 

where 



is the Hamiltonian. 

Similarly, in the flavor state basis, one has 

^|^/W=^/^/W, (12) 

where 

In the mass eigenstate basis, the matter term (the potential matrix) is 

Vm = U-WfU. (13) 

Thus, the total Hamiltonian in the mass eigenstate basis is given by 

J^m = Hm + U'^VfU. (14) 

Equation (|ll]) has the solution 

Mt) = e-'^"Vm(0). (15) 

Inserting t = L into Eq. (0), one finds the solution in the mass eigenstate basis to be 

ijm{L) = e-'^-'^^miO) = UUL)tpUO) (16) 

and in the flavor state basis 



^ f//(L)^;(0). (17) 



Inserting the potential matrix in the flavor basis and the CKM matrix into Eq. (|T3D 
yields 

/ f/^i UelUe2 UelUeS \ 
V^ = A\ Ue2Uel U^2 Ue2Ue3 , (18) 

\Ue3Uel f/e3f/e2 f/fs / 

which is a (real) symmetric matrix. Then, the Hamiltonian in the mass eigenstate basis is 

/Ei+AU^, AUelUe2 AUelUeS \ 

M'„, = AU,2Uel E2 + AUl^ AUe2Ue3 . (19) 

V AU,sU,, AU,sU,2 Es + AU^J 

Of course, we also have Uf{L) = e~*'^/^, but we prefer to first work out Um{L) = e"*-^^-^, 
since the structure of J^m is much simpler than that of Jiff, and afterwards perform the 
transformation to Uf{L). 

Next, we need to find the exponential of the matrix —iJifmL in order to obtain Um{L). 
Generally, when one wants to find the exponential of an A^ x A^ matrix M, one has to 
calculate 



00 _ 
^" = E>"- (20) 



The characteristic equation of the matrix M is given by 

X(A) = det(M - XI) = X^ + cn-iX^'^ + ■ ■ ■ + ciA + cq = 0, (21) 

where Cj,z = 0,l,...,A^ — 1, are some coefficients. Cayley-Hamilton's theorem implies that 
A in the characteristic equation can be replaced with the matrix M itself, giving 

M^ + cn-iM^~^ + ■■■ + ciM + Co/ = 0, (22) 

where / is the N x N identity matrix. Hence, one has 

M^ = -CN-iM^'^ ciM - Co/, (23) 

which means that 

MP = c%l^M^-^ + ■■■ + cf'^M + cfl, for any p>N, (24) 

where cf , 2 = 0, 1,... ,A^ — 1, are some coefficients. The exponential of the matrix M can 
then be written as 

N-l 

e^ = aol + aiM + ■■■ + qn-iM^'^ = J2 ^"^"' (^5) 

n=0 



where ai, i = 0,1, ... ,N — 1, are some coefficients. Thus, the infinite series in Eq. ( pOD is 
reduced to only A^ terms, where N is the dimension of the matrix M. 
An arbitrary N x N matrix M can always be written as 

M = Mo + ^(trM)/, (26) 

where Mq is an A^ x A^ traceless matrix, i.e., trMo = 0. Splitting up M as in Eq. (^) is 
useful, since the identity matrix / is commuting with all other matrices. 

Now using Eqs. p5|) and (p6|), one can write the exponential of the matrix —iJ^mL as 

g-i^™L = (f)e-'^^ = (f) [aol + ai{-iLT) + a2{-iLTf] = {aol - iLTai - L'^T'^ai) , (27) 

where = e"*^*''-^™/^ is a complex phase factor, T = M'm — (tr M'm)l j'i is a traceless matrix, 
and Oq, a\, and a^ are coefficients to be determined. The infinite series in Eq. (pO]) is here 
reduced to only three terms. The trace of the matrix M'm is 

tr ^^ = El + ^2 + ^3 + A. (28) 

The matrix T can now be written as 



r = (r, 



ah) 



A[/2, - iA + i (Ei2 + Eis) ^f/elf/e2 AV ,^U , 



AU,^U,2 AUl^-\A+\{E2i + E2^) AU,2U, 



AU,,U,s AU,2Ue3 At/23 - |A + 1(^31 + ^32, 

(29) 

where Eat = Ea — Ef,. The six quantities Eat, where a, 6 = 1, 2, 3 and a ^ b, are not linearly 
independent, since they obey the following relations: 

Eba = —Eab, (30) 

Eu + E23 + E31 = 0. (31) 

This means that only two of the EabS are linearly independent, e.g., E21 and £"32. 

The following linear system of equations will determine the coefficients oq, ai, and a2 
introduced in Eq. (^Tf ) above, 

g-iLAi =aQ- iLXiai - L^Xla2 

e-*^^2 = ao - iLX2ai - L'^\la2 , (32) 

g-iLA3 =aQ- iLXsai - L'^Xla2 

where Xa, a = 1,2, 3, are the eigenvalues of the matrix T, i.e., solutions to the characteristic 
equation 

A^ + C2X^ + ciA + Co = 

with 

C2 = -trT = 0, 



Cl — ^11X22 — T12T21 + TiiT; 



To'iT: 



23 J 32 



M' 



33 ~" ^13^31 + ^22^33 - ^^- .,^ 
\A [Ul (E21 + E31) + f/e2 (^12 + ^32) + Ul (Ei3 + E. 



~ 9 (-^12 + -^13 + -^23 + E12E13 + E21E23 + £'31-^32) , 



23 



Co = - det T = - Ayl^ + 1^2 [^yfi (i?21 + i?3l) + f/e2 (^12 + i?32) + f/'s (^13 + ^^23 

— g^ft/g]^ {E21 + E23) (-E'31 + -£^32) + U^2 (-^12 + -E'13) (-E'31 + E32) 

+ ^e3 (-^12 + -^13) (-£'21 + -E'23) + I {E12 + -^13 + -E23 + E12E13 + E21E23 + E31E32) ] 

— 27 (-^12 + -^13) (-^21 + -E'23) (-E'31 + -E'32) • 



Note that the coefficients cq, ci, and C2 are all real. Introducing the relativistic limit (see 
Appendix 0), E21 and ii^32 can be written in terms of the mass squared differences Am^ 



Am^J and AM^ 



Am^2) ^s 



E. 



21 



Aw? 
~2E 



and -£"■ 



32 



AM2 



where i? is the neutrino energy. 

The solutions to the characteristic equation are |l^ 



Ai = -- (si + S2) - -C2 + i— (si - S2) = -- (si + S2) + %— (si - S2) , (33) 



"2 (si + S2) - 2C2 - «— (si - S2) 



1 a/S 

'2 (si + S2)-«^(si-S2), (34) 



A3 = Si + S2 - -C2 = Si + S2, 



(35) 



where 



and 



„2^ 1/21 1/3 



,, = [r + (g3 + r^)!/^] ^'^ , S2 = [r - {^ + r^)^/^] 



2\l/2ll/3 






3-i 9-2 3^1? 
'^ = 6('^1'^2 ~ 3Co) — 27C2 = ^2*^0" 

Since T is a Hermitian matrix (T is a real symmetric matrix), it holds that g'^ + r^ < and 
all eigenvalues Aq, where a = 1, 2, 3, are real. Furthermore, inserting q = ci/3 and r = — co/2 
into the inequality g'^ + r^ < 0, one finds cf /27 + Co/4 < from which one can conclude that 
Cl < and |ci| > 3|cop''^/v^, since Cq,Ci E M. Thus, the roots (eigenvalues) can be written 
as 



Ai = -A - 



--C1C0S 



1 / 1 

- arctan — 

3 \ Co 




+ V^ 



sm 



1 / 1 

- arctan — 
3 \co 



-cl - ^c? 



(36) 



1 



-Ci COS 



1 / 1 

- arctan — 
3 \ Co 



27 



V— Ci sin 



1 / 1 

- arctan — 

3 \ Co 




(37) 



A, 



-ci cos 




(38) 



Observe that Ai + A2 + A3 = 0. This relation is satisfied because of the fact that the trace 
of T is zero. In addition, the relations A1A2 + A1A3 + A2A3 = ci < and A1A2A3 = — cq are 
fulfilled. 

The Aa's are, up to a constant, the energy eigenvalues of the neutrinos in the presence 
of matter. The measurable quantities are the energy differences |Aa — Af,|, which in the 
relativistic limit are related to the effective matter mass squared differences Am^j, as 



An — Ah 



2E 



a, b 



1,2,3, 



a ^ b, 



(39) 



where E is again the neutrino energy. 

Introducing the following matrix form representation for the system of equations in 
Eq. (HD, 



Aa, 



(40) 



where 



-iLXi 
-iL\2 



one obtains the solution 



A 




1 


-iLXi 


-L^A? 


1 


—iL\2 


-L^A^ 


1 


-iLXs 


-L'Xl 



and 



tL^ {XIX3 - AiA^) tL^ {XiXl - XIX2 




-L' (Ai - XI) 
-iL (Ai - A3) 




A^A3) 



+ e 



-JLA2 



(A?A3 



'^^' (Ai - Ai) + e-^^^2 (Ai - A?) + e 
'"^^^ (A3 - A2) + e-^^^2 (Ai - A3) + e 



L' (A? - XI)' ) f e 
2L(A2-Ai) / Ve 

-AiAi) + e-*^^«(AiA2-A?A2 
'^^^ (A2 - Ai 



(41) 



where 



D = det A = iL^ (AiA^ - X\X2 + A2Ai - AiA3 + A3A? - A^Ai) . 
Thus, the exponential of the matrix —iJ^^L in Eq. (^) can be written as 



^ ^ D 



■j [e~^^^^ (A^A3 - A,A^) + e-^^^ (A^A^ - A^As) 

+ 6-^^"== (A?A2 - AiA^) ] / + [e-^^^i (A^ - A^) + e"^^^^ (A? - A^) 
+ 6-'^^== (A^ -Xl)]T+ [e-^^^^ (A3 - A2) + e-'^^^ (A3 - Ai) 

+ e-'^^^(Ai-A2)]T4. (42) 

Inserting the expression for D into Eq. (021) gives 



e-'=^™^ = TT ^4^ ^0e-^^^i [A2A3/ - (A2 + A3)T + T 

[Ai — A2)[Ai — A3) 

^ (A3-A.)(A3-A.) ^'"'"' ['''^' - ''' ^ '^'^ ^ ^^^ • <'" 

Using the various relations for the eigenvalues finally yields 

,-^™l _ .-.LX-. (A^ + ci)/ + AiT + T^ , -..A. (Aj + cQJ + A2T + T^ 
^ 3A2 + C1 ^^ 3Ai + ci 

^^ 3Ai + ci ' ^ ^ 

which can be written as 

^ 1 

tt 3A2 + ci L 

The evolution operator for the neutrinos in the flavor basis is thus given by 

3 1 . . 

f//(L) =e-*^^^ = ^7e-^^'"^?7-^ =0^e-*^^"— (A^ + ci)/ + A^t + t^ , (46) 



a=l 



3X1 + ci 



where T = UTU^^ . Equation ( ^Uf ) is our final expression for Uf{L). 

Let us pause for a moment to contemplate Eqs. ( ^5]) and (|46|). Since .^ = f/,^^^~^, 
it is clear that T = J^f — (tr=^)//3 = J^j — (tr ,^^)J/3 due to the invariance of the trace 
under transformation of U. In fact, the characteristic equation is also invariant under U 
and therefore so are the coefficients cq, ci, C2, and the eigenvalues Ai, A2, A3. However, the 
expression for Jiff is much more complicated than that for J^m, which is the reason why we 
work with Jifm instead of J^j. 

The formula (^) is our main result for the evolution operator. It expresses the time (or 
L) evolution directly in terms of the mass squared differences and the vacuum mixing angles 
without introducing any auxiliary matter mixing angles. 
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IV. PROBABILITY AMPLITUDES AND TRANSITION PROBABILITIES 

The probability amplitude is defined as 

A^p = {l3\Uf{L)\a), a,f3 = e,n,T. (47) 

Inserting Eq. ( ^BD into Eq. (^) gives 

A., = <Pl^e ^^^^^ , (48) 

a=l " 

where 5aj3 is Kronecker's delta and the entries Taj3 and (T'^)ai3, a,P = e,fi,T, are given in 
Appendix g Here f^p = fp^ and {f'^)af3 = {T'^)iia- 

The probability of transition from a neutrino fiavor a to a neutrino fiavor (3 is thus 
defined by the expression 

Pap = \Aai3t = Al^f^A^p, a,P = e,n,T. (49) 

Inserting the expression for the probability amplitude into the definition of the transition 
probability, one finds 



a=l 6=1 



-iLjXa^X,,) (-^g + Cl)^a/3 + A^Ta/j + (T )a/3 (A^ + Ci)^a/3 + AfeTg,/; + (T )q,/3 

3A2 + ci 3X1 + ci 



(50) 



Since Sap, Tap, and (T^)ai3 all are symmetric, it holds that Pa/3 = Ppa- 
Setting A = gives 

{h\U^{L)\a) = (6|e-^^™^|a) = e-'^"%b, a, 6= 1,2,3, (51) 

3 
Aap = {l3\Uf{L)\a) = Y,UaaUpae-'''^\ a,f3 = e,^,T, (52) 

a=l 

and 

3 3 

Paf3 = Sa/3-4:^ ^UaaUf3aUabUf3b siu^ Xab, a, P = 6, fl, T, (53) 

a=l 6=1 
a<b 

where Xab = EabL/2, which are the well-known transition probabilities for vacuum oscilla- 
tions (see, e.g., Ref. [Q). 

We can write the probabilities for oscillations in matter in a form analogous to the ones 
for the vacuum probabilities, as 

D / V^ l/^a "I" '^lJ"a/3 "I" ^a-'a/S "I" l-i Ja/3 

^«/3 - 2^ 



' jXl + C,)6af3 + Xgfap + jf^ ' 

3A^ + C, 



_ . ST^ \^ (Aa + Cl)5a/3 + XaTa/3 + (T )a/3 (A^ + Ci)5q,/3 + XbTa/3 + (T )a/3 . 2 ~ 

a=l 6=1 '^^'^ + ^1 "^^b + ^1 

a<6 

a, /3 = e, /i, r, (54) 

11 



where Xab = (Aa — \b)L/2. For L = Eq. ( ^61) reads 

I = E TV2^l(^l + ci)I + Xaf + n (55) 

,=1 3A, + ci 

This means that the corresponding matrix element version of the above formula is given by 



3 



Sap=> ^^T^^ -, a,(3 = e,^i,T. (56) 

a=l ^^- + ^1 

Hence, the transition probabilities in matter are 

D X A \^ \^ (-^a "'" '^l)'^Q/3 + XaTap + (T )„/3 (A^ + Ci)(5q,/3 + AftTa^ + (T )q,/3 . 3 ~ 

p., - 5., - 4 ^ ^ ^3^^-^^^ ^^^^^ sm X.,, 

a,P = e, n,T. (57) 

From unitarity, there are only three independent transition probabilities, since the other 
three can be obtained from them, i.e., from the equations 

Pee + Pefi + Per = 1, (58) 

P^e + ^MA. + Pf^r = 1, (59) 

P.e + Pr^. + Prr = 1, (6O) 

where Pe^ = P/^e, Per = Pre, and P^r = Pr;^. We will choose Pe^, Per, and Pfj,r as the three 
independent ones. The expression for the transition probabilities in Eq. ( ^Tj) thus simplifies 
to 



^«/^ = -4 2^ 2^ ,X2^., ,X2^.. sm x,„ (61) 

a=l b=l 
a<b 



3A2 + ci 3A2 + ci 



which is valid only for a ^ (3. 

V. DETERMINING THE AUXILIARY MATTER MIXING ANGLES 

The flavor states can be expressed as linear combinations of either the vacuum mass 
eigenstates {A = 0) in the basis Hm or the matter mass eigenstates {A ^ 0) in the basis 
l-Lu- The corresponding components are related to each other as follows: 

^f = U^Pm, (62) 

V^/ = f/^'^M, (63) 

where U = U{6i, 62,03) and U^ = 1/^(6^,62^,6^'^) are the unitary mixing matrices for 
vacuum and matter, respectively. Here 6f^ , i = 1,2,3, are the auxiliary (effective) matter 
mixing angles. 
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Combining these expressions for the flavor components, one obtains the following relation 
between the two different sets of mass eigenstate components: 



IpM = Rlp-n 



(64) 



where 

R = {U^'Y^U. 

The matrix R is, of course, a unitary matrix (even orthogonal, since U and U^ are real). 
This means that the matter mixing matrix can be expressed in the vacuum mixing matrix 

as 



f/*^«,^f ,^3^0 = u{e,,e2,e:,)R'\ 



(65) 



The relations between the different bases can be depicted as in the following commutative 
diagram: 



^m e Hr. 



R 



i^M e n 



M 



u 



u 



M 



Hf^^j 



From this diagram one readily obtains the Hamiltonian in presence of matter J^^ 



M as 



which is diagonal in the mass eigenstate basis in matter, TIm- 
Due to the invariance of the trace, we have 



Tm = M'm - |(tr JTm)/ = RTR-^, 



where Ti 



M 



is a traceless diagonal matrix. This implies that 



-iLXa 



a=l 



6A^ + Ci 

which is actually equivalent to the system of equations (|32D. 
Inverting the first relation above, we obtain 

Uf{L) = (l)Ue-'^^U'^ = (l)UR-^e-'^^''RU-^ = (j)U^^e-'^'^"{U^'^)~\ 



(66) 



(67) 



(68) 
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From this one immediately obtains the transition amphtude 

3 

A^p = {(3\Uf{L)\a) = (Pj^UaaU^ae-'''^^ a,(] = e,fi,T, (69) 

a=l 

where the f/^'s are the entries in the matrix U^ . This allows us to identify the U'^ matrix 
elements with those in Eq. (H3) as 



jMtjM ^ (-^g + Cl)5ap + XgTa p + (T j^p 

"" ^'^ 3A2 + ci 



'Jca'Jl3a = 5T^- , a,fJ = e,^,T, = 1,2,3. (70) 



The matter mixing angles can therefore be expressed as follows 



ef = arctan ^ = arctan ^3^-^^ + (^')eA^ /^^^ 

U^ A3re. + (T2)e.' ^ ^ 



' Ai + Ci + A3Tee + ( T2^ 

3Ai + ci 

. + {f^)ea 3Af + Ci 



6*2 =arcsm(7g3 =arcsm^/ ^^^ , (72) 

?3^ = arctan -^-^ = arctan '^ _'^" ' ^ , ^^° " '^ ' "^ , a = fj.,T. (73) 



t/if ^ A2Te, + (T2),„3A? + Ci 



The transition probabilities in matter can thus be written as 

3 3 

PaP = S^p-4Y,Y.^aaUZU^bU^b^in'Xat, «, /? = e, /X, T. (74) 

a=l 6=1 
a<b 

The above formulas for the transition probabilities are thus in effect exactly the same as the 
formulas in Eq. (pTj). 

When the matter density goes to zero, i.e., A — > 0, then U^^ — > Uaa (6'f^ -^ Gi) and 
Xab — > Xab- Also Eqs. ([57| ) and ( |7iD will become identical to Eq. ( |5^ ) when A = 0. 

To connect this treatment to the usual treatment in terms of the auxiliary matter mixing 
angles, we have to find the matrix R. Using Eq. (pBD, the matrix i?, which diagonalizes the 
matrix T, can be constructed from the eigenvectors of the matrix T as 



R=\^l\. (75) 



I.e., 




i?-l = i?r = ( vi V2 V3 ) , (76) 
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where 



Here 




1,2,3. 



a„ 



Ti2 {T33 — Xa) — TiaT^ 



23 



T23 [Tu — Xa) — T12T] 
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Pa 



Tl2 (^33 '~ ^a) — Ti^T'^ 



23 



Tl3 (T22 — Xa) — T12T:; 



23 



V^^l + W+^'^ 



1,2,3. 



The orthogonal matrix R is of course diagonahzing both T and J^m- 
Equation (^) can now be written in matrix form as 



u^l uM uM 

TjM TjM TjM 

'^/il '^/i2 '^/iS 

jjM jjM jjM 

^t1 ^t2 '-^t3 



MnM 



X 



C^'C. 



•^3 ^2 



5f 



■^MnM qM QMnM nMnM qM qM qM cMr^M 

'-'1 ^2 ^3 "-"l ^2 



On O1 0-\ On ^3 ^1 ^3 

qM qM cMr^MryM cMr'M oM cM/^M r^M r^M 

'-'1 '-'3 ~ '-'2 ^1 ^3 ^'-'1 ^3 ~ '-'2 '-'3 ^1 ^1 ^2 

C'2'-^3 *J3'-^2 ^2 

— 03O1 — 010203 O1U3 — 01O203 0102 
01O3 — 02O1O3 — 01O3 — 02O3O1 O1O2 

aiZi 02-22 as-zs 
(5izi P2Z2 P3Z3 

Zl Z2 Z3 

Uel Ue2 Ue3 \ ( (Viji (V2J1 [y3h 
Uf,l U^2 U^3 (Vi)2 (V2)2 (V3)2 

Url Ur2 Ur3 J \ (vi)3 (V2)3 (V3)3 




(77) 



where S^^ = sin 6^^ and C^^ = cos6f^ for i = 1,2, 3. We can then obtain the matter mixing 
angles from the above matrix relation as 



Of 



M 



arctan 



Sf 



M 



CI 



M 



arctan 



(-5'3Ci - 5'i5'2C3)Q^3^3 + (C1C3 - 5'i5'25'3)/33Z3 + S1C2Z3 
{S1S3 — S2CiC3)a3Z3 + { — S1C3 — S2S3Ci)(33Z3 + C1C2Z3' 



(78) 



nM ■ oM 

62 = arcsm 63 



arcsin(C2C3a32;3 + S3C2P3Z3 + S2Z3), 



(79) 



lAf 



arctan 



SI 



arctan 



^2^-30^2^2 + *-'3C2P2^2 + 02^2 

C2C3aiZi + S3C2I31Z1 + 5*2^1 ' 



(80) 



The conditions for resonance are given by the extremal points of the probability ampli- 
tudes. These will coincide with extremal points for the modulus of the amplitudes, thus for 
the products U^lJj¥ or for the expressions in the right-hand side of Eq. ( [TOD as functions of 
A. 
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VI. DISCUSSION 

The main result of our analysis is given by the time evolution operator for the neutrinos 
when passing through matter with constant matter density in Eq. p6|) expressed as a finite 
sum of elementary functions in the matrix elements of the Hamiltonian. Also the transition 
probabilities in Eq. (|6T|) and the neutrino energies in presence of matter Ea = Xa, where 



a = 1,2, 3, given by Eqs. (pBD -fpSD belong to our main results as do the expressions for 
the matter mixing angles in Eqs. (|7T])-(^). In our treatment these auxiliary matter mixing 
angles play no independent role and are not really needed. It is convenient, though, to 
express the matter mixing angles in the matrix elements of the Hamiltonian to look for the 
resonance conditions. 

As an illustration of the resonance phenomena, we have plotted the energy differences 
|Aa — Afe|, where a,b = 1,2,3, a ^ b, and the behavior of the matter mixing angles 9f^ , 6*^^, 
9^ as a function of the matter density A for two different cases. First for two mass squared 
differences that are close and bimaximal mixing [^, i.e., two large mixing angles and one 
small. The results are presented in Figs. |I] and ^. Then for two widely separated mass 
squared differences and all the mixing angles small. These results are presented in Figs. |^ 
and Ij. The energy differences \Xa — \b\ are related to the effective matter mass squared 
differences Am^^ by the following relations 

\Am1ij\ = 2E\\a — \b\, a, 6=1,2,3, a ^ b. 

As can be seen from Figs. |I|-^, the resonances occur (i.e., the mixing is maximal) when 
the energy levels in presence of matter approach each other as function of the matter density 
A. As the mass squared differences approach each other in magnitude, the resonances move 
closer together and the mixing angle 6f^ , which looks quite flat (constant) at large separations 
of the mass squared differences, will have a modest increase above the location of the (higher) 
resonance of the mixing angle 62^ . The smaller the vacuum mixing angles are, the shaper are 
the resonance peaks. For bimaximal mixing the resonance occurs only for the small mixing 
angle, but even for matter densities above the resonance region the values of the auxiliary 
mixing angles are changed appreciably from their vacuum values, as is shown in Fig. ^. 

Our results for the auxiliary matter mixing angles are in agreement with the previous 



approximative and numerical calculations in Refs. [P,|^,|TT],|T^,[T^ . 

We have also calculated the transition probabilities Pap for neutrino oscillations in matter 
as functions of the matter density A. The results are illustrated in Figs. ^ and |^. These two 
figures correspond to the only physically measurable quantities. Our illustration is chosen 
to correspond to values of the mass squared differences that are close to those obtained 
from analyses of atmospheric neutrino data. In Fig. |^(c), we can see a sharp drop in the 
probability Pgg at a matter density of about 5- 10^^"^ eV (same parameter values as in Figs. |I| 
and 1^) for a sensitivity of L/E ^ 1.6 ■ 10^ eV"^ as the neutrinos pass through the Earth. 
This corresponds to an enhancement of Pe/j, and Per at the same value of A, whereas the 
probability P^,- wiggles in the same region. Due to unitarity only three of these probabilities 
are of course independent. 

The Cayley-Hamilton formalism could in principle be used to derive neutrino oscillation 
formulas in matter for arbitrary numbers of neutrino flavors. However, the calculations for 
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four or even more flavors will be much more tedious than in the case of three neutrino flavors. 
For four or more flavors, the unitary mixing matrices have more complicated structures than 
in the three flavor case. Another complication, for more than four flavors, is that one will 
not, in general, be able to find analytical solutions to the characteristic equation, since only 
characteristic equations up to degree four are analytically solvable. 

In Appendix 0, we have given, as a comparison and reference, the derivation of the neu- 
trino oscillation formulas in matter for two neutrino flavors also using the Cayley-Hamilton 
formalism. We observe that the derivation for two flavors is much easier to make than that 
for three flavors. The evolution operator Um{L) can of course be written in a much more 
compact form for two flavors than for three; compare Eqs. ( ^Sj) and (|C5D . 
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APPENDIX A: INTRODUCING THE RELATIVISTIC LIMIT 

We will here introduce the relativistic limit by reparametrizing the expressions for the 
energy differences. Any two of E21, -E32, and E13 can be chosen as independent variables. 
Thus, we use Eab = Am^^/(£'a + Eb) for a,b = 1, 2, 3 and a ^ b, where Am^^ = m^ — m^. 
Without lack of generality, we will assume that there is a mass ordering with m-i < 1712 < m^. 
Apart from the masses, we will take the momentum as a common independent variable and 
assume that the momentum is the same for all components of the mixed state. However, 
it is convenient to use instead the energy as an independent quantity. This choice can 
again be made in several ways. For definiteness, we will here choose as physical variables 
E={Ei + £'2)/2, Am^ = Am|i, and AM^ = Am|2. This gives 

Aw? , , , 

S. = ^, (Al) 




!S.vn?^ 



Since in the applications I\MjE ^ 1 (^ 10^^), we find the excellent approximation 

E^2 ^ ^^. (A3) 

Equations ( |Al| ) and (|A3|) lead to the following expressions for the coefficients C2, Ci, and Cq: 

C2 = 0, (A4) 
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ci ^ -\a^ + -4 [Ul {AM' + 2Am') + UJ^ {AM' - Am') - f/fg {2AM' + Am')] 
^ (AM^ + Am^ + AM' Am') , (A5) 

-^A' + ^ [f/fi (AM^ + 2Am') + UJ, {AM' - Am') - U^ {2AM' + Am')] 

A 

— — [U^^ {2AM' + Am') {AM' - Am') + UJ^ {2AM' + Am') {AM' + 2 Am') 
Ul^ {AM' + 2Am') {AM' - Am') - {AM^ + Am^ + AM' Am') ] 
—^ (2AM' + Am') (AM' + 2Am') (AM' - Am') . 



+ 



(A6) 



APPENDIX B: THE ENTRIES OF THE MATRICES f AND f^ 

The entries of the 3x3 real symmetric matrix T' are 

{T')n = i [A' (f/fi + l)+2A {Ul, - 1) {E,, + E,,) + \ (Eu + Ers)'] , (Bl) 

{T')22 = I [A' {Uj^ + \)+2A {Ul^ - \) (^21 + E,^) + 1 (^21 + E,;)'] , (B2) 

(T^)33 = I [A' {Ul, + \)+2A {Ul, - \) (^31 + E,,) + 1 (^31 + E,^'] , (B3) 

(T2)i2 = {T')2l = lUelUe2A (A + E^ + E^^) , (B4) 

(T2)i3 = (T2)31 = |f/elf/e3^ {A + E12 + ^32) , (B5) 

(T2)23 = {T')^2 = \Ue2Ue^A (A + E21 + ^3l) • (B6) 

For the matrices T and T^ we have the following expressions for the entries: 

3 3 

Tk/3 = ^ ^ UaaUiSbTab, a, P = 6, fl, T, (B7) 



and 



a=l 6=1 



3 3 



(f2)„^ = ^^t/,,t/^,(T2),,, «,/? = e,/x,r. (B8) 



a=l fe=l 

For a ^ [3 Eq. ( p7| ) simplifies to 



3 



^Q/3 = 2_^ UaaUfiaEa- (B9) 



a=l 



APPENDIX C: THE TWO FLAVOR NEUTRINO CASE 

Considering only a two flavor neutrino oscillation model instead of a three neutrino, the 
free Hamiltonian in the mass eigenstate basis and the potential matrix are 

^1 ^ _, ,^ f A 

'2 



H^=[ Q^ and V>=( J. (CI) 



In the mass eigenstate basis, the potential matrix is given by 

Vm = U-'VfU, (C2) 

where 

rr = f ^^^ ^'"'^^ = f ^°*^ ^ sin 6* \ 
[U^i U^2j ~ \-sin9 cosOj 

is the usual 2x2 orthogonal matrix, which describes neutrino mixing with two neutrino 
flavors. Thus, we have 

and 

^ _ TT ,T/ _ f El+ AUl^ AUelUe2 \ /p.N 

-^m-^m+Vm- y AUelUe2 ^2 + AU^^ ) ' ^^ > 

Again, we have to find the evolution operator Um{L) = e"*'^™^. Similar but simpler formulas 
than in the three flavor neutrino case now give 



' V e-^^^" -^ {\J + T)=<P f cos(LA) - il sm{L\)T] 

a=i ^^'^ V A y 



(C5) 
where = e-'^^''^^/'^, T = .JTm - (tr JC)//2, Ai = -A, A2 = A, and A = V- det T. Here 



A 



1.2 ^21/, ^. a;2i\ _ E21L Am^L ^ 2_ 2 2 

-A^ - — {^Acos2e- —J, 3:21 = ^— = ^^ , Am = m2 - m^. 



The quantity A comes from the solutions, Ai and A2, to the characteristic equation A^ + 
ciA + Co = 0, where ci = -trT = and cq = detT = -^74 + AU^^E2i/2 + AUI2E12/2 + 
E12S21/4 = -AV4 + 0:21 {A cos 2^ - a:2i/L) /L. 

The evolution operator Um{L) can now be written in matrix form as 

f/™(^) = 0(^ f*), (C6) 

where a = cos(LA) - i (Acos2^ - 2x2i/L) sin(LA)/(2A) and (3 = -iAsin2^sin(LA)/(2A). 
Similarly, the evolution operator Uf{L) is 
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Uf{L) = UU^{L)U-' 

cos 9 sinO \ f a P \ l' cos 9 — sin ^ 
— sin 61 cos 6 J \ (3 a* J \ sin 6 cos 6 

a cos^ 9 + a* sin^ 9 + psm 29 ^^ sin 29 + p cos 29 
^^ sin 29 + (3 cos 2^ a sin^ 9 + a* cos^ ^ - /5 sin 20 

The e/i-element of the Uf{L) matrix is 



(C7) 



(y — (y Toi 

(f//(L))e^ = sin20 + /5cos20 = -i0^sin20sin(LA). (C8) 

Taking the absolute value and squaring Eq. (|C^), we obtain the transition probability Pe^ 
as 

Pe, = I (f^/(^))e,r = -^^sin'29sin\L~X) 



(Lxy 

sin2(LA) 



X21 sin^ 29 



L2 [A74 - X21 {A cos 20 - xai/L) /L] 

sin^(LA) 



1 + AL [AL /{AX21) - cos 20] /a;2i 
sin^ 20 



sin^(LA), (C9) 



sin^ 20 + [cos 20 - AL/{2x2i)y 
where LX = \/ x\^ — ALx2\ [cos 20 — y4L/(4x2i)]. Introducing the definitions 

sin^ 20^ ^ — ""^^^ . (CIO) 

sin2 20 + (cos20-A/E2i)^ 

and 



X2x = \\x\- ALx2i(cos29-- — j = X2i'i/sin2 20 + | cos20 - — j , (Cll) 

we can write 

Pe^ = sin2 20^sin2x2i. (C12) 

One could also immediately identify the definitions ( |C10| ) and ( |C11| ) from Eq. ( |C8| ) by 
observing that the factor 0:21 sin 20/ (LA) should correspond to sin 20^^ and that the argument 
LA in the factor sin(LA) should correspond to X2\- This identification corresponds to the 
one in Eq. (|70D above in Sec. 0. The factor — i0 is just a complex phase factor. 
Setting 

cos20=^ (C13) 

-E21 
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in Eqs. (|CTg) and (|CT|) gives 

sin^ 29 = 1 (maximal mixing) 



(C14) 



Pel, = sin 



3^21 + 



AL 



AL 



a;2i - 



sin X21 W 1 



A 



E. 
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A 



E. 



21 



(C15) 



where Eq. ( pi3| ) is the so-called MSW resonance condition. 

Inserting A = into Eq. ( P12|) leads, of course, back to Pg/i 
transition probability formula for z/e-z/^ oscillations in vacuum. 



sin^ 2^sin^ X21, the old 
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FIG. 1. The differences |Aa — A^l, a,b = 1, 2, 3, a ^ b, as a function of the matter density A for 
(9i = 6*3 = 45° (bimaximal mixing), 02 = 5°, Am^ = 10"^ eV^ AM^ = lO"-'^ eV^ and E = W GeV. 
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FIG. 2. The quantities sin^ 261, 
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1,2,3, as a function of the matter density A for 



0^ = 03 = 45° (bimaximal mixing), 02 = 5°, Am^ = 10"^ eV^, AM^ = IQ-^ eV^, and ^ = 10 GeV. 



Note that sin 20^ = 1 corresponds to maximal mixing. 
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FIG. 3. The differences |Aa — Af,|, a,b = 1, 2, 3, a ^ b, as a function of the matter density A for 
= 6)2 = 6*3 = 10°, Am2 = 10-10 eV^, AM^ = IQ-^ eV^, and E = W GeV. 
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FIG. 4. The quantities sin"' 26'-" , i 
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sin2 2^F 



1,2,3, as a function of the matter density A for 
6*3 = 10°, Am2 = 10-1° eV^, AM^ = IQ-"^ eV^, and ^ = 10 GeV. Note that 
1 corresponds to maximal mixing. 
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FIG. 5. The transition probability Pge ^& a function of the matter density A for 
= 6*3 = 45° (bimaximal mixing), 0^ = 5°, Am^ = lO""^ eV^, AM^ = lO'^ eV^, and 



LjE = r]{2R(Q/W GeV) ~ 3.23 x 10"^ eV"^ rj, where i?® ~ 6378 km ~ 1.62 x 10^^ eV- 
the (equatorial) radius of the Earth, (a) rj = 1/5, (b) r/ = 1, and (c) r] = 5. 
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FIG. 6. Transition probabilities as functions of the matter density A, using the same parameter 
values as for Fig. |(c). (a) Pee, (b) Pe^, (c) Per, and (d) P^r- 
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